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ABSTRACT 


Polarization  of  light  due  to  metallic  reflection  is  at  best  only  partially  explained  by  current 
models.  A  new  model  is  proposed  to  account  for  this  phenomenon  involving  "native  clusters". 
These  clusters  are  assumed  to  be  an  integral  part  of  metallic  bulk  surfaces  and  to  manifest  their 
dielectric  properties  in  the  visible  region.  An  index  of  refraction  may  be  computed  based  on  a 
pseudo-Brewster  angle,  that  angle  of  incidence  where  the  minimum  perpendicular  reflectance 
occurs.  Next  the  model  identifies  the  cluster-to-bulk  surface  area  ratio  of  the  material  and 


formulates  the  contribution  of  reflection  from  these  areas  by  utilizing  the  Fresnel  equations. 
Reflectance  from  the  remainder,  the  bulk-surface  area,  is  computed  by  applying  conventional 
metallic  reflection  techniques.  By  combining  the  two  contributions  of  reflectance  for  angles  of 
incidence  from  0  to  90  degrees,  reflectance  curves  are  generated  for  gold,  silver,  nickel  and 
aluminum.  These  calculations  assume  a  wide-band  radiating  source  of  light.  Although  the 
magnitude  of  reflectance  for  the  model  showed  some  variance  over  the  range  of  angles  of 
incidence  (within  10%),  the  curves  are  similar  in  shape.  The  minimum  perpendicular  reflectance 
and  grazing  angles  of  incidence  are  consistent  with  experimental  findings.  The  computed  surface 
ratio  for  each  metal  varies  indirectly  with  the  metals’  normal  reflectance  value,  as  expected  from 


the  model.  Spectral  information  on  the  complex  index  of  refraction  is  included  to  infer  probable 


cluster  size. 
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I.  DESCRIPTION  OF  MECHANISMS  FOR  METALLIC  REFLECTION 


A.  INTRODUCTION 

Optical  polarization  was  discovered  by  E.  L.  Malus  in  1808.  However,  it  was  not 
until  the  nineteenth  century  and  the  formulation  of  electromagnetic  theory  developed  by 
Maxwell  that  this  phenomenon  could  be  described  analytically,  at  least  partially. 
Moreover,  the  theory  addresses  the  results  more  than  the  mechanism.  J.  Bennett  and  E. 
Bennett  from  the  Naval  Weapons  Center  are  quoted  in  Ref.  1  as  saying  in  reference  to 
the  electromagnetic  theory, 

"This  theory  is  phenomonological  in  that  instead  of  trying  to  explain  why  materials 
have  certain  fundamental  characteristics,  it  concentrates  on  the  resulting  properties 
which  any  material  with  those  characteristics  will  display.  Electromagnetic  theory 
has  little  or  nothing  to  say  about  why  a  material  should  have  these  particular  optical 
properties  ..." 

This  is  a  strong  statement  concerning  the  present  lack  of  understanding  involved  in  the 
polarization  process  during  metallic  reflection.  Bennett  and  Bennett  continue  to  state  that 
at  the  present  time  solid  state  physics  can  only  partially  explain  the  mechanism  which 
causes  polarization  for  some  surfaces.  In  this  chapter,  currently  accepted  mechanisms 
are  mentioned  and  the  concept  of  "native  clusters"  will  be  introduced. 

B.  PRESENT  SCOPE  OF  RESEARCH 

Clusters  are  a  collection  of  atoms  which  have  properties  between  bulk  and 
molecular  entities.  Clusters  are  also  known  as  microparticles,  or  finite  systems.  They 
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are  primarily  studied  in  the  free  state  [Ref.  2J.  "Native  clusters"  [Ref.  3],  refer  to  the 
covalent  aggregation  of  a  fixed  number  of  atoms  at  the  surface  energy  level  of  the  metal. 
This  clustering  is  presumed  to  exist  at  all  times,  although  the  cluster  site  position  on  the 
surface,  varies  with  time. 

Having  performed  a  thorough  literature  search,  and  having  been  sponsored  by  the 
Naval  Postgraduate  school  to  attend  an  international  symposium  on  the  physics  and 
chemistry  of  finite  systems  [Ref  4],  the  author  has  found  that  presently  there  is  intense 
interest  within  the  scientific  community  concerning  clusters.  However,  despite  the  wide- 
ranging  topic  material  presented  at  the  conference  [Ref.  4],  there  is  a  noticeable  lack  of 
research  involving  bound  "native  clusters"  and  their  impact  on  surface  science.  It 
appears  that  the  main  thrust  at  present,  is  on  defining,  building  and  manipulating  the  free 
clusters,  which  are  created  through  vaporization  of  bulk  material  by  high  energy  lasers. 
Since,  as  stated  in  Refs.  2  &  4,  free  clusters  are  derived  from  the  bulk  material  and  by 
themselves  display  properties  other  than  the  bulk  from  which  they  were  produced  and 
that  many  clusters  assembled  together  begin  to  take  on  the  properties  of  the  bulk,  it 
seems  logical  that  research  involving  "native  clusters"  would  add  to  the  understanding 
of  how  clusters,  a  distinct  phase  of  matter  [Ref  2],  interact  within  metallic  surfaces  and 
modify  their  emission  behavior. 

C.  MECHANISMS  OF  POLARIZATION 

Currently  there  are  two  accepted  models  which  partially  describe  the  behavior  of 
reflected  waves  from  metallic  surfaces:  The  free  electron  theory  which  is  best  suited  for 
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the  infrared  spectrum  and  below  and,  conversely  the  quantum  model  works  well  in  the 
ultraviolet  region  of  the  energy  spectrum  and  beyond. 

1.  Free  Electron  Theory 

The  free  electron  theory  is  based  on  a  model  where  actually,  two  types  of 
electrons  reside  on  the  surface.  Bound  electrons  or  electrons  which  are  covalently 
bonded  to  the  atoms  which  make  up  the  crystalline  surface  of  the  metal  and  free  electrons 
which  are  associated  with  good  conductors  and  are  free  to  move  about  in  a  random 
pattern  thereby  producing  no  net  current  [Ref.  5].  Introducing  the  complex  index  of 
refraction,  the  interaction  of  the  free  electrons  with  an  electromagnetic  wave  can  be 
accounted  for  by  the  imaginary  part,  which  in  turn  depends  on  the  properties  of  the 
metal.  More  specifically,  on  the  conductivity  as  it  relates  to  the  number  of  free  electrons 
in  the  metal  [Ref.  5].  However,  as  one  will  note  in  Chapter  III,  the  model  is  limited  to 
the  portion  of  the  spectrum  in  the  infrared  (IR).  However,  two  important  facts  are 
validated  by  this  model. 

•  Electromagnetic  waves  in  the  near  infrared  and  longer  wavelengths  are  exclusively 
absorbed  at  the  surface  by  interaction  with  the  free  electrons,  (see  Chapter  III) 

•  Magnitudes  of  parallel  and  perpendicular  polarized  reflectance  can  be  obtained 
from  the  Fresnel  equations  over  the  same  portion  of  the  energy  spectrum  given 
above. 

2.  Quantum  Mechanics 

For  reasons  explained  in  Chapter  III  the  quantum  model  is  best  suited  for 
explaining  the  mechanism  in  the  ultraviolet  (UV)  and  higher  energy  portion  of  the 
spectrum.  High  energy  waves  possess  enough  energy  to  influence  conduction  bands,  it 
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is  theorized,  above  the  Fermi  level  but  below  the  work  function  level.  If  this  is  the  case 
then  it  may  be  possible  to  associate  these  bands  with  dipoles  which  are  the  necessary 
ingredient  for  polarization  to  occur.  However,  this  linkage  is  not  well  developed  at  this 
time  due  to  the  complexity  of  the  interaction  between  high  frequency  waves  and  the  band 
electrons  [Ref.  5]. 

3.  Polarization  Due  To  "Native  Clusters" 

"Native  clusters"  are  postulated  to  exist  as  discrete  sites  on  the  surface,  to  be 
transparent  and  non-absorbing  in  the  visible  portion  of  the  spectrum  [Ref.  6].  The  theory 
that  "native  clusters"  do  exist  on  metallic  surfaces  is  proposed  in  Ref.  3.  In  this  work, 
Bibiarz  complements  the  free  electron  model  used  in  describing  thermionic  emission  from 
metals.  It  is  postulated  [Ref.  3],  that  clusters  are  intrinsic  regions  which  exist  on  the 
surface  of  bulk  metal  and  which  alter  all  emission  properties  of  the  bulk  (i.e.,  electron, 
atom  and  photonic).  A  thermionic  model  is  developed  which  includes  "native  clusters" 
as  a  percentage  of  the  total  surface  area  of  the  metal.  Combining  the  bulk  contribution 
to  the  contribution  from  the  clusters,  the  total  thermionic  emission  is  determined  over  a 
range  of  temperatures.  As  shown  in  Ref.  3,  this  model  better  predicts  thermionic 
emission  than  the  pure  free  electron  model  in  the  case  of  refractory  metals. 

The  conjectured  shape  of  clusters  is  interesting  from  the  standpoint  of  their 
integration  into  the  bulk  surface.  Clusters  are  mostly  surface  area.  In  fact  Ref.  2  states 
that  tightly  packed  clusters  of  20  atoms  may  have  only  one  atom  in  their  interior.  This 
would  lead  one  to  think  of  clusters  as  having  much  more  surface  than  depth  when 
composing  a  portion  of  the  bulk  metal.  Further  hampering  the  defining  of  the  cluster 
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shape  is  the  inability  to  observe  cluster  atomic  structure  directly  with  current  electron 
microscopy  methods  if  the  cluster  is  less  than  about  two  nanometers  in  size  [Ref.  4], 
However,  some  cluster  shapes  are  being  inferred  from  mass  spectral  data  and  magic 
numbers  associated  with  the  material  [Ref.  4]. 

If  the  cluster  size  is  around  one  micrometer  across  [Ref.  3],  then  one  can  discount 
from  being  applicable  to  this  model  those  portions  of  the  spectrum  which  are  made  up 
of  wavelengths  either  too  big  to  affect  the  cluster  or  are  so  small  that  their  interaction 
with  the  electrons  in  the  cluster  is  negligible.  For  these  portions  of  the  spectrum,  either 
the  free  electron  model  for  long  wavelengths  or  the  quantum  model  for  short  wavelengths 
may  be  a  more  satisfactory  explanation  of  the  mechanism. 

Free  metallic  clusters  can  be  created  through  vaporization  of  the  bulk  metal  by  high 
energy  lasers.  Clusters  created  in  this  manner,  when  studied,  do  not  have  the  same 
optical  properties  as  the  bulk  [Ref.  2].  Therefore,  it  is  reasonable  to  assume  that  clusters 
residing  as  a  portion  of  the  total  metallic  surface  would  not  share  the  host  optical 
properties  and,  in  fact,  could  have  the  properties  stated  at  the  beginning  of  this  chapter. 
Hence,  it  follows  that  if  the  clusters  do  exhibit  dielectric  behavior,  their  presence  can 
be  accounted  for  in  the  theory  of  electromagnetic  radiation  and  by  the  Fresnel  equations. 
By  adding  this  portion  to  the  contribution  from  the  bulk  (free  electrons),  total  metallic 
reflectance  may  thus  be  determined.  As  such,  a  mechanism  whereby  waves  of  the  near 
IR  and  visible  portion  of  the  sp>ectrum  are  polarized,  is  established. 
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n.  BACKGROUND  ON  ELECTROMAGNETIC  WAVE  BEHAVIOR 


When  an  electromagnetic  wave  encounters  a  boundary  between  two  different  optical 
media,  such  as  an  air-dielectric  system,  it  is  well  known  that  the  wave  is  resolved  into 
two  component  parts,  one  being  the  reflected  wave  the  other  being  the  transmitted  wave 
[Ref.  7].  In  this  work  we  consider  "white"  light  and  unpolarized  electromagnetic  waves, 
and  as  such,  several  properties  of  the  wave  as  it  encounters  a  dielectric  boundary  can  be 
explained  by  electromagnetic  theory  based  on  Maxwell’s  equations  [Ref.  5]. 


A.  INDEX  OF  REFRACTION 

The  speed  of  light  in  vacuum  is  defined  and  given  by  Equation  2.1  below. 


c  = 


1 


(2.1) 


Where  pto  is  the  permeability  of  free  space  and  Eq  is  the  permitivity  of  free  space. 

For  optical  materials,  the  index  of  refraction  is  defined  to  be  the  ratio  between  the 
speed  of  light  (c)  in  a  vacuum  and  the  speed  of  light  in  the  medium  (v)  [Ref.  7]  and  is 
given  below  in  Equation  2.2. 


n  = 


c 

V 


(2.2) 


However,  for  the  derivations  that  follow  it  will  be  more  correct  to  make  a  substitution 
for  V  in  Equation  2.2  as  follows:  It  is  understood  [Ref.  5]  that  for  good  dielectrics  the 
magnetic  permeability  {fi)  of  the  dielectric  is  equal  to  the  magnetic  permeability  of  free 
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space.  Therefore  it  follows  from  electromagnetic  theory  that  the  velocity  of  the  wave 
(light)  in  the  dielectric  as  given  by  Equation  2.3  [Ref.  7]. 

(2.3) 

Finally,  by  substitution  of  Equations  2. 1  and  2.3  into  Equation  2.2  the  index  of  refraction 
for  dielectrics  is  given  by  Equation  2.4. 


For  the  remainder  of  this  paper  the  permitivity  term  (^)  will  refer  to  a  relative  permitivity 
defined  as  the  ratio  of  permitivity  of  the  dielectric  over  permitivity  of  free  space.lt  should 
be  noted  that  the  index  of  refraction  for  dielectrics  is  a  real  number  unlike  the  index  of 
refraction  for  conductors  which  is  given  in  complex  form.  The  reason  for  this  difference 
will  become  clear  later. 

The  index  of  refraction  is  useful  in  defining  a  relationship  between  the  angle  of 
incidence  (<f))  of  the  incident  wave  and  the  refraction  angle  (0')  of  the  refracted  wave 
(Figure  2.1  )  known  as  Snell’s  law  [Ref.  7]. 


sin<))  _ 
sin<J)'  n 
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Figure  2.1  Diagram  Of 
Reflection  And  Refraction 
Between  Two  Optical  Media  n 
and  n* 


B.  FRESNEL  EQUATIONS  FOR  DIELECTRICS 

The  description  of  the  Fresnel  equations  for  the  transmitted  and  reflected  light 
waves  is  straightforward  and  summarized  below.  A  detailed  derivation  can  be  found  in 
any  book  on  electromagnetic  wave  propagation,  such  as  Ref.  5. 

Once  a  light  wave  encounters  a  dielectric  boundary,  the  wave  is  split  into  two  new 
waves.  From  Figure  2. 1 ,  one  wilt  note  there  is  a  reflected  wave  and  a  transmitted  wave 
each  with  a  parallel  and  perpendicular  polarization  component  associated  with  it.  If  we 
let  the  amplitudes  of  the  incident,  reflected,  and  transmitted  waves  be  denoted  by  I,  R, 
T,  respectively,  as  shown  in  Figure  2.1,  then  apply  the  boundary  condition  that  the 
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tangential  components  be  continuous  across  the  boundary  [Ref.  8],  the  results  are  given 
by  Equation  2.6  [Ref  5]. 

Since  the  assumption  was  made  in  Ref.  8  that  the  clusters  behave  as  a  dielectric,  n'  = 
J  e.  Making  this  substitution  into  the  above  equation  and  solving  for  Rj^,  R|,  T^,  T|  the 
well  known  Fresnel  equations  for  reflectance  are  obtained  and  given  below  by  Equation 
2.7. 


COS<j)'  T^  =  cos4>(A,  -  i?,) 

ye7  T^  =  +  i?,) 

r,  =  A,  +  R, 

y^cos(j)'  r,  =  y/tl  (.fl,  +  R,) 


RSDP 

RSDS 


2 

tan(4)i  -  (j)')  f 

tan  (4)_j  +  V) 

2 

sin((J)_^  -  <1)0 

A, 

sin  {({)_;  +  <1)0 

(2.6) 


(2.7) 


Note  that  RSDP  and  RSDS  are  reflectance  coefficients  and  are  the  result  of  squaring  the 
magnitudes  of  the  intensities  of  the  reflected  and  incident  waves,  as  in  evaluation  of  the 
pointing  vector,  then  dividing  by  the  square  of  the  component  amplitude  of  the  incident 
wave  which  normalizes  the  equations. 


C.  BREWSTER’S  LAW  FOR  DIELECTRICS 

Brewster’s  law  simply  states  that  there  exists  a  specific  angle  of  incidence  which 
is  material  dependent,  such  that  the  reflected  light  from  the  dielectric  is  plane  polarized 
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(i.e.,  RSDP  =  0).  From  the  top  formula  in  Equations  2.7  one  will  note  that  this  can 
only  occur  when  the  angle  between  the  reflected  and  transmitted  waves  is  90°  (<f),  +<j)’ 
=  90°)  and  hence  the  tangent  goes  to  infinity.  The  definition  of  the  Brewster  angle,  leads 
to  Equation  2.8.  This  angle  will  be  relevant  to  the  development  of  the  "native  cluster" 
model,  in  Chapter  IV. 

tanO  =  —  (2.8) 

n 

D.  REFLECTION  AND  REFRACTION  PROPERTIES  OF  A  CONDUCTOR 

When  an  electromagnetic  wave  impinges  on  a  metallic  boundary  (i.e.,  a  good 
conductor)  a  surface  current  is  immediately  established  .  Unlike  dielectrics,  whose 
surface  structure  is  composed  of  bound  atoms  and  electrons,  a  conductor’s  surface  is 
usually  composed  of  small  non  orientated  crystals  containing  freg  electrons.  These  free 
electrons  will  move  so  as  to  cause  a  surface  current  when  in  the  presence  of  an  external 
electric  field. 

1.  Complex  Index  Of  Refraction 

According  to  [Ref.  5],  by  replacing  the  real  index  of  refraction  (n)  as  given 
in  Equation  2.4,  to  a  complex  index  of  refraction  (N)  defined  by  Equation  2.9  an 
elementary  model  can  be  developed  whereby  the  optical  constants  e,  a,  fx,  are  utilized 
as  in  Equation  2.10,  it  defines  the  real  and  imaginary  parts  of  the  complex  index  of 
refraction.  Note  that  a  is  the  conductivity  of  the  metal. 
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n  -  ik 


(2.9) 


N  = 


n 


2  _ 


2  .2  t  4  11^0^  ^ 


e 


n^k'^  = 


2  ,2  t  4  |i"g2  _ 


e 


(A  e 


H  e 


(2.10) 


Work  published  in  1903  by  E.  Haygen  and  H.  Rubens  [Ref.  5]  indicates  that  the  results 
of  Equations  2.10  are  accurate  for  wavelengths  of  10  cm  (infrared  waves)  or  longer. 
This  however,  will  lead  to  some  error  between  the  calculated  and  experimental  results 
of  these  optical  properties  in  the  visible  to  ultraviolet  portion  of  the  spectrum. 
Therefore,  this  model  needs  to  be  restricted  to  near  infrared  wavelengths. 
Notwithstanding  the  aforementioned,  the  complex  index  of  refraction  arises  when 
dielectric  electromagnetic  theory  is  expanded  to  encompass  conductors.  Figure  2.2  from 
Ref.  1,  shows  the  trends  in  the  reflectance  curves  for  increasing  n  and  k.  It  is 
particularly  helpful  in  showing  that  wave  propagation  into  the  conductor  is  small,  and 
hence  can  be  neglected.  By  defining  the  "skin-depth"  (d)  as  the  distance  an 
electromagnetic  wave  will  travel  into  a  conductor  before  dissipation  of  it’s  energy  to  a 
value  of  1/e,  and  linking  this  to  the  inverse  of  the  extinction  coefficient  [Ref.  5],  one  can 
show  that  the  distance  of  wave  propagation  into  the  medium  is  given  by  Equation  2.11. 


d 


c 

v/8npou) 


(2.11) 
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Figure  2.2 


Reflectance  curves  for  various  values  of  n  and 
k  derived  from  the  free  electron  model. 


12 


Angle  of  incidence,  deg 


From  Equation  2.11,  one  will  note  that  high  conductivity  (on  the  order  of  10'^  sec 
Gaussian  units)  and  high  frequencies,  as  in  visible  waves,  will  result  in  a  very  small  skin 
depth.  As  proof  an  example  of  how  little  the  wave  will  be  transmitted  into  the 
conductor.  Ref.  5  gives  a  skin  depth  of  d  =  10'^  cm  for  infrared  wave  into  copper  a 
typical  metallic  conductor.  This  fact  plays  an  important  role  in  the  development  of  the 
bulk  reflection  model  as  outlined  in  Chapter  IV,  since  we  can  consider  the  metals  to  be 
opaque  in  the  visible  portion  of  the  spectrum. 
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ra.  FREE  ELECTRON  AND  QUANTUM  MODELS 


A.  THE  FREE  ELECTRON  MODEL  OF  METALLIC  REFLECTION 

In  Chapter  II,  wave  reflectance  from  a  conductor  was  addressed  using  the  classical 
electromagnetic  theory  based  on  Maxwell’s  equations,  together  with  the  introduction  of 
the  complex  index  of  refraction.  In  this  section  the  free  electron  model  (proposed  by  P. 
Drude)  will  be  shown  and  comments  on  its  short  comings  in  dealing  with  the  higher 
frequency  wave  reflectance  will  be  made. 

The  basis  of  the  theory  for  the  free  electron  model  is  that  in  a  good  conductor  two 
types  of  electrons  must  reside  at  the  surface.  The  first  type  are  called  bound  electrons, 
or  those  which  are  held  in  a  quasi-rigid  state  by  the  parent  atom,  similar  to  the  dielectric. 
The  second  type,  called  free  electrons,  is  a  well  established  model  for  conductors,  in  that 
the  surface  is  full  of  these  electrons  which  are  not  bound  to  any  atom  and  can  move 
fairly  freely  about.  As  discussed  in  the  previous  chapter,  the  complex  index  of  refraction 
is  frequency  dependent.  However,  Ref.  5  states  that  the  complex  conductivity  and 
emissivity  are  both  frequency  dependent,  which  is  shown  by  Equations  3.4  and  3.6. 

1.  Derivation  Of  The  Complex  Conductivity  Term 

From  statistical  methods  of  the  theory  of  kinetics  of  gases,  the  motion  of  the 
free  electrons  can  be  described  by  a  model  of  a  proportional  damping  force  acting  in  the 
opposite  direction  of  the  free  electron  motion  [Ref  5].  This  motion  which  is  cyclic  in 
nature  due  to  the  electric  field  being  cyclic,  can  be  described  by  the  standard  second 
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order  differential  equation  of  motion  for  a  damped  spring  mass  system  under  the 
influence  of  the  forcing  function  E  as  shown  below. 


mr  +  /7?Pr  =  eE  (3.1) 

The  solution  of  Equation  3.1  is  given  in  complex  form  below  (i  =  imaginary  index). 

f  =  - - ® -  E  (3.2) 

/n(o)2  +  ipo)) 

Where  r  is  the  position  vector  of  the  periodic  motion  of  frequency  (co)  the  electron  in  an 
electric  field  with  a  damping  coefficient  (/3).  The  decay  time  (t)  is  inversely  related  to 
the  damping  coefficient  (/3)  which  is  on  the  order  of  lO'''*  sec  [Ref  5].  Differentiating 
Equation  3.2,  to  obtain  the  velocity  of  the  electron,  and  substituting  this  into  Equation 
3.3  below 


2. 


Complex  Emissivity 


since 


As  stated  in  Chapter  II  for  conductors,  the  index  of  refraction  is  complex  and 


n 


pe' 

N 


(3.5) 


It  follows  that  if  the  index  of  refraction  is  complex  then,  by  definition  as  shown  in 
Equation  3.5,  so  must  the  emissivity  be  complex.  For  low  frequencies  Ref.  5  states  that 
the  real  part  of  £  is  small  compared  to  the  imaginary  part  and  hence  Equation  2.9 
becomes  the  expression  for  the  complex  £. 


g  2 


1  - 


/no)  ( (i)  -  i  P ) 


(3.6) 


Hence  the  dependence  on  frequency  is  again  established  as  it  was  for  the  conductivity. 

3.  Analysis  Of  The  Free  Electron  Model 

The  free  electron  model  as  presented  predicts  well  the  optical  constants  of  the 
conductor  for  long  wavelengths  (IR,  near  or  far),  but  does  not  accurately  predict  the 
behavior  of  the  constants  and  hence,  the  reflection  of  the  wave,  in  the  visible  portion  of 
the  energy  spectrum.  Evidence  of  this  is  presented  in  Ref.  5.  The  problem  stems  from 
the  fact  that  Equation  3.6  as  stated  above  is  derived  based  on  low  frequency  electric 
fields,  in  order  to  eliminate  the  contribution  of  the  bound  electrons.  However,  the 
contribution  from  the  bound  electrons  becomes  more  significant  as  the  frequency 
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increases  thereby  leading  to  the  breakdown  in  the  model.  It  is  stated  in  Ref.  5  that  the 
free  electron  model  can  not  be  modified  in  order  to  extend  it’s  usefulness  into  the  higher 
frequency  ranges  due  to  the  complexity  of  the  interaction  between  the  electric  field  and 
the  metal  in  these  ranges.  It  refers  the  reader  to  the  theory  of  quantum  mechanics  for 
an  explanation. 

B.  QUANTUM  THEORY 

Band  theory  of  solids  is  presented  in  various  references  (see,  for  example,  Ref  9). 
For  metals,  the  theory  seems  most  appropriate  in  the  violet  and  smaller  wavelengths. 
Interband  electron  transitions  in  the  near  infrared  are  best  known  in  semiconductors.  The 
quantum  mechanical  expression  for  the  polarizability  is  of  the  same  form  as  that  obtained 
from  the  classical  (Drude)  model. 

Weak  interband  transitions  appear  to  decrease  the  reflectance  of  aluminum  around 
1.4  eV  .  In  silver,  the  reflectance,  well  into  the  "vacuum"  ultraviolet  region,  is 
characteristic  of  interband  transitions  spread  over  a  wide  spectral  region.  At  3.9  eV  a 
sharp  decrease  in  reflectance  identifies  the  free-electron  plasma  frequency.  Plasma 
oscillations  are  a  non-quantum,  collective  oscillation  of  the  free  electrons.  Discussion 
of  the  present  state  of  this  field  is  beyond  the  scope  of  this  work. 
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IV.  MODEL  DEVELOPMENT 


A.  "NATIVE  CLUSTER"  DESCRIPTION 

"Native  clusters"  are  postulated  to  be  clusters  which  form  not  as  a  result  of  the 
presence  of  "foreign"  atoms  deposited  on  the  surface,  but  rather  because  of  an 
equilibrium  of  surface  atom  configurations  [Ref.  3].  The  "native  cluster"  model  (Figure 
4.1)  is  able  to  account  for  a  number  of  the  "nonconductor"  properties  of  reflected  light 
from  metallic  surfaces  since  in  a  cluster,  the  electrons  remain  attached  to  the  atoms 
which  make  up  the  cluster  and  hence  do  not  contribute  to  the  Fermi  level.  Biblarz 
postulates  in  Ref.  6  that  a  patchwork  of  "native  clusters"  that  exist  on  virgin  metallic 
surfaces  (based  on  work  published  in  Ref.  3)  does  alter  polarization  upon  metallic 
reflection  .  The  dispersement  of  such  "native  clusters"  over  the  bulk  surface  may 
produce  dispersed  dipoles  thereby  providing  the  surface  with  "nonconductor"  attributes 
as  exhibited  by  polarization,  by  the  reflection  of  visible  light  from  the  surfaces,  as  well 
as  other  attributes. 

B.  FORMULATION  OF  EQUATIONS 

Based  on  Ref.  3,  there  is  assumed  to  exist  an  equilibrium  surface  area  ratio  of 
cluster  to  bulk  material  unique  to  any  given  metal.  Therefore  reflected  light  will  have 
contributions  from  the  bulk  (rg)  and  clusters  (rg).  The  total  surface  area  (Sj)  is  defined 
as  the  sum  of  the  cluster  area  (Sg)  and  the  bulk  area  (Sg)  as  shown  in  Equation  4.1. 
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(4.1) 


Ss  Sq  Sj 

Dividing  Equation  4. 1  by  the  total  surface  area  (normalizing)  and  defining  the  fraction 
Sj/St  as  X,  the  cluster  to  total  surface  ratio,  one  can  substitute  this  value  into  Equation 
4.1  and  solve  for  the  normalized  bulk  to  total  surface  area  ratio  as  in  Equation  4.2. 


(1  -  X) 


(4.2) 


Total  reflectance  (t-^)  is  then  given  in  [Ref.  5]  as 

Tj.  =  (1  -  xlr^  +  xr^  (4.3) 

Since  the  bulk  can  be  considered  opaque  in  the  near  IR  and  visible  portion  of  the 
spectrum  [Ref.  6]  (the  skin  depth  being  on  the  order  of  10’  cm),  the  contribution  from 
the  bulk  is  simply  the  reflectance  of  the  bulk  times  it’s  normalized  surface  area. 
However,  clusters  behaving  as  dielectrics,  are  presumed  to  be  transparent  and  will  allow 
a  portion  of  the  waves  to  be  transmitted  through  the  cluster  and  then  reflected  by  the 
metallic  substrate  while  the  remainder  is  reflected  directly  at  the  surface.  As  such, 
reflectance  from  the  clusters  maybe  composed  of  direct  reflectance  (rsd)  and  indirect 
reflectance  (rjj)  leading  to 
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A  gray  body  is  defined  to  have  characteristics  of  a  black  body  where  the 
absorptance  is  equal  to  the  emittance.  Assuming  that  the  metallic  surface  behaves  as  a 
gray  body  (a  good  assumption  for  infrared  to  visible  wavelengths)  then  the  absorptance 
must  equal  the  emissivity  (Kirchhoff  s  law)  [Ref.  1]  as  shown  in  Equation  4.5. 

a(k,T)  =  e(X,T)  (4.5) 

Energy  flux  incident  on  surface  can  be  reflected,  transmitted,  or  absorbed.  Applying 
Equation  4.5  to  the  flux  balance  equation  [Ref.  1]  the  relationship  between  absorptance 
{a  or  e),  reflectance  (r)  and  transmittance  (t)  is  shown  to  be 

e{X,^,T)  +  r(A.,<t),T)  +T(A.,(j),r)  =  l  (4.6) 

C.  BULK  REFLECTANCE 

Transmittance  of  the  light  wave  into  the  bulk  is  shown  to  be  negligible  as  stated  in 
Ref.  5,  hence  Equation  4.6  may  be  written  as 

Tg  =  1  -  e  (A.,  r)  cos(|)  (4.7) 

The  angle  of  incidence  (<f»)  is  defined  as  the  angle  between  the  incident  ray  and  the 
normal  to  the  surface  (see  Figure  2.1). 

According  to  Ref.  6  the  emittance  is  well  represented  in  the  IR  region  of  the  spectrum 
by  Equation  4.8, 


e(X,D  -  36. 


(4.8) 


[onm-cm)  / A.  (\iw) 

where  p  (a  function  of  temperature)  is  the  resistivity  of  the  metal  in  ohm-cm  and  X  is  the 
wavelength  of  the  incident  wave  in  pm.  The  non-spectral  form  of  the  emmittance  is 
determined  by  integrating  over  the  entire  spectrum  and  shown  in  Equation  4.9. 

, _  (4.9) 

£  (T)  =  0.51 52 yjp  (ohm-cm)  -T^iK) 

D.  CLUSTER  REFLECTANCE 

Native  clusters  possessing  dielectric  properties  will  cause  incident  waves  to  either 
be  reflected  directly  (rs<i)  or  indirectly  (rgd  after  transmission  through  the  cluster  and 
reflection  from  the  metallic  substrate  (see  Figure  4.2).  As  shown  by  Equation  4.4. 


1.  Direct  Cluster  Reflectance 

Direct  cluster  reflectance  is  calculated  using  the  Fresnel  equations  for 
dielectrics  [Ref.  7].  The  reflected  light  has  two  components,  parallel  and  perpendicular, 
representing  the  polarization  effect  upon  reflection  from  a  dielectric. 
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2.  Indirect  Cluster  Reflectance 


The  indirect  portion  of  the  reflection  is  much  more  difficult  to  define  as  it 
requires  knowledge  of  the  shape  and  size  of  clusters.  Some  insight  regarding  the  size 
of  the  "native  cluster"  is  discussed  in  Chapter  VI.  No  attempt  was  made  to  include  this 
contribution  but  preliminary  estimates  indicate  that  the  magnitude  of  the  contribution  to 
reflection  will  be  on  the  order  of  20%  or  less. 
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V.  COMPUTER  MODEL  DESIGN 


A.  COMPUTER  ALGORITHM  DEVELOPMENT 

Utilizing  equations  derived  in  Chapter  IV,  a  computer  algorithm  named 
"REFLECTA"  (see  Appendix  A)  was  developed  to  calculate  the  parallel  (RTOTLS)  and 
perpendicular  (RTOTALP)  reflectance  from  the  metallic  surface  for  angles  of  incidence 
from  0°  to  90°.  Reflectance  curves  were  generated  from  the  data  by  invoking  "Easy- 
Plot",  a  locally  developed  software  package  used  with  FORTRAN  VII  at  the  Naval 
Postgraduate  School,  and  plotting  RTOTLP  or  RTOTLS  against  the  angle  of  incidence 
((t>)  for  each  metal. 

B.  CURVE  HTTING 

White  light  reflectance  curves  for  gold,  and  silver  based  on  experimental  results 
[Ref.  7],  were  considered  to  be  the  baseline  from  which  the  "native  cluster"  model  could 
be  compared  to,  then  scrutinized  for  similarities  and  differences.  Being  able  to  produce 
reflectance  curves  from  the  model  which  matched  the  experimental  curves,  would  be 
proof  for  the  validation  of  the  model. 

Three  parameters  ,  bulk-cluster  index  of  refraction  (n'),  cluster  to  bulk  surface  area 
ratio  (x)  and  emittance  (t)  were  utilized  in  obtaining  the  best  curve  fit  for  a  particular 
metal. 
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C.  EQUATIONS  USED  IN  COMPUTER  ANALYSIS 


Initially  the  parameters  were  fixed,  based  on  results  reported  in  Ref.  7  for  gold 
and  silver,  which  measured  reflectance  from  a  "white"  tungsten  light  bulb,  probably 
operating  at  3000  K.  Plots  for  these  curves  are  shown  in  Chapter  VI  as  Figures  6. 1  and 
6.2.  However,  curves  were  also  generated  for  aluminum  and  nickel  in  order  to  expand 
the  data  base  for  future  trend  analysis  (see  Figures  6.3  and  6.4).  Referring  to  the 
REFLECTA  program,  four  terms  were  considered  a  critical  part  in  plotting  of  the 
reflectance  coefficient  curves,  and  are  defined  below  for  continuity  purposes  as  they  will 
be  repeatedly  referred  to  in  the  remainder  of  the  thesis. 

•  RTOTLP  is  the  total  surface  parallel  reflectance  coefficient  of  the  metallic  surface. 

•  RTOTLS  is  the  total  surface  perpendicular  reflectance  coefficient  of  the  metallic 
surface. 

•  RSDP  is  the  parallel  reflectance  coefficient  contribution  from  direct  cluster 
reflection  (Equation  5.1). 

•  RSDS  is  the  perpendicular  reflectance  coefficient  contribution  from  direct  cluster 
reflection  (Equation  5.2). 

RSDS  and  RSDP  were  calculated  by  Equations  5.1  and  5.2. 


RSDS 


sin  (4)  -  4)0 
sin{4»  4>0 


(5.1) 


RSDP 


tan  (4)  -  4>0 
tan  (4>  -*■  4>^) 


(5.2) 


From  Equations  5.3  and  5.4,  RTOTLP  and  RTOTLS  were  calculated. 
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RTOTLS 


(1  -  X)RB  +  xRSDS 


(5.3) 


RTOTLP  =  (1  -  x)RB  *  xRSDP  (5.4) 

It  was  noted  that  the  angle  of  incidence  where  RTOTLP  was  a  minimum  occurred 
about  four  degrees  right  of  the  experimental  curve  given  in  Ref.  7.  By  noting  that  this 
minimum  is  associated  with  the  pseudo-Brewster  angle  for  the  metal,  as  defined  by  the 
minimum  perpendicular  reflectance  as  shown  in  Ref.  7,  calculation  of  the  index  of 
refraction  for  the  bulk-cluster  aggregate  (n')  is  determined  by  the  following  equation. 

n'  =  (5.5) 

where  <f>^  is  the  pseudo-Brewster  angle  and  n,  is  the  index  of  refraction  of  air  or  vacuum 
(equal  to  1.0).  However,  by  adjusting  the  incident  angle  to  match  the  experimental 
pseudo-Brewster  angle  the  entire  reflectance  curve  was  shifted  up  or  down  depending  on 
the  direction  of  shift  in  </>  caused  by  the  matching  of  the  pseudo-  Brewster  angle.  At  this 
point  it  became  clear  that  a  sensitivity  analysis  for  each  parameter  was  required  in  order 
to  formulate  a  systematic  approach  whereby,  the  three  parameters  could  be  varied  in  a 
convergent  manner  thus  providing  the  best  curve  fit  by  the  model. 

D.  SENSITIVITY  ANALYSIS  OF  PARAMETERS 

As  stated  above  in  Section  B  of  this  chapter,  due  to  the  interdependence  of  the 
three  parameters  used  in  the  generation  of  the  curves  and  their  complexity  it  is  necessary 
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to  perform  an  analysis  on  each  parameter  in  order  to  determine  what  effect  a  small 
change  in  each  has  on  the  reflectance  curves. 

1 .  Procedure 

Sensitivity  testing  was  performed  using  gold,  by  varying  one  parameter  in 
small  increments  while  holding  the  other  two  constant.  Total  reflectance  was  recorded 
at  two  angles  of  incidence  each  chosen  for  their  critical  nature  in  the  curve.  The  first 
was  at  </>  =  2°  where  the  incident  light  is  virtually  normal  and  hence  the  cos((/>)  is  nearly 
one  thereby  eliminating  the  dependence  of  angle  of  incidence  at  that  point.  The  second 
is  located  at  a  point  where  RTOTLP  is  at  a  minimum  since,  this  will  provide  the  best 
information  on  curve  shifts  due  to  small  parameter  changes  at  the  Brewster  angle.  It  is 
felt  that  the  best  curve  fit  would  be  obtained  if  these  points  could  be  made  to  match  with 
experimental  results. 

2.  Findings  From  Sensitivity  Analysis 

The  results  of  the  sensitivity  testing  are  shown  in  Figure  5.1  and  are 
summarized  below. 

•  Small  changes  in  the  parameters  produced  linear  increments  of  change  in  the 
curves. 

•  Increasing  x  decreased  RTOTLP  at  0  =  2°  and  the  minimum  RTOTLP.  There 
was  no  change  in  which  was  expected  since  this  is  governed  by  n’  only. 

•  Varying  e  had  no  effect  on  the  shape  of  the  curves.  Thus  this  parameter  can  be 
determined  by  Equation  4.9  and  fixed. 

•  Increasing  n’  increases  RTOTLP  at  <^i  =  2°;  produces  no  change  in  RTOTLP^; 
and  increases 
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3.  Partial  Conclusions  From  Sensitivity  Analysis 

From  the  displayed  linearity  above  and  the  equation  of  a  line  we  can  generate 
two  simultaneous  equations.  Letting  y  be  the  value  of  RTOTLP  at  an  angle  of  incidence 
equal  to  2°,  we  can  solve  for  the  two  unknowns,  slope  of  the  line  (m)  and  the  intercept 
(b)  for  a  given  n',  set  initially  by  the  pseudo-Brewster  angle.  Once  this  is  accomplished, 
the  value  for  x  can  be  determined  for  that  metal  by  setting  y  equal  to  the  experimental 
value  of  RTOTLP  at  normal  incidence  and  solving  for  x.  The  curves  are  then  generated 
using  these  parameter  values.  The  angle  of  incidence  corresponding  to  RTOTLP„i„  is 
compared  with  the  experimental  value  and  a  small  change  is  made  in  n'  in  order  to  shift 
the  curve  towards  the  experimental.  This  new  n'  is  then  used  in  determining  a  new  value 
for  X  and  the  process  repeats  itself  until  <l>am  experimental  matches  with  that  of  the 
program.  It’s  important  to  note  at  this  time  that  the  initial  estimate  of  x  should  be  chosen 
so  as  to  comply  with  current  wisdom  regarding  the  area  ratio  of  cluster  to  bulk. 
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VI.  RESULTS 


A.  CURVE  nX  EXPERIMENTAL  VERSUS  THE  MODEL 

Utilizing  the  aforementioned  "native  cluster"  model  and  the  computer  algorithm 
developed  in  the  previous  chapter,  total  reflectance  curves  were  generated  for  four 
metals,  gold,  aluminum,  nickel,  and  silver  (Figures  6.1  thru  6.4).  These  data  as  well 
the  results  produced  by  the  model  are  for  wide-band  light  (i.e.,  white  light)  where  the 
source  temperature  was  assumed  to  be  3000  K.  Experimental  values  for  gold  and  silver 
are  shown  as  x’s.  The  definition  of  these  parameters  is  necessary  in  order  to 
incorporate  the  experimental  curves  displayed  in  Ref.  7,  for  comparison  purposes  to 
those  of  the  model. 

Comparison  of  the  reflectance  curves  produced  by  the  model  to  the  experimental 
curves  given  in  Ref.  6  reveal  the  following  similarities  listed  below. 

•  At  normal  incidence  (4>  =  0”),  both  parallel  and  perpendicular  reflectance  have 
the  same  magnitudes.  Further,  this  common  magnitude  is  in  agreement  with  the 
value  measured  by  experimental  means. 

•  Perpendicular  reflectance  increases  as  the  angle  of  incidence  increases 
throughout  the  range  of  values  in  a  similar  fashion  to  the  experimental  values. 

•  Parallel  reflectance  decreases  to  a  minimum  located  at  an  angle  of  incidence 
associated  with  the  pseudo  Brewster  angle  of  the  metal  as  does  the  experimental 
value  of  parallel  reflectance. 

•  Parallel  reflectance  increases  for  angles  of  incidence  greater  than  the  pseudo- 
Brewster  angle  to  the  grazing  angle  of  incidence  {<)>  =  90°). 
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TOTAL  REFLECTANCE  (GOLD) 


Figure  6.1  Reflectance  Curves  Of  Gold,  x's  correspond  to 


neasurements  from  Ref.  7. 
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TOTAL  REFLECTANCE  (SILV 


Figure  6.2  Reflectance  Curves  Of  Silver.  x's  correspond  to 
laeasureraents  from  Ref.  7. 
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TOTAL  REFLECTANCE  (NICKEL) 
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Figure  6.3  Reflectance  Curves  Of  Nickel. 
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TOTAL  REFLECTANCE  (ALUMINUM) 
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•  At  the  grazing  angle  of  incidence,  both  parallel  and  perpendicular  reflectance  are 
equal  to  unity  as  are  the  respective  experimental  values. 

The  cluster  model  is  intended  to  permit  agreement  between  experimental  and  cluster 
model  values  of  reflectance,  for  two  angles  of  incidence,  (j>  ranging  from  0°  to  90°,  as 
well  as  the  angle  of  incidence  where  the  perpendicular  reflectance  is  at  a  minimum. 
However,  for  intermediate  angles  of  incidence  the  reflectance  values  are  determined  by 
the  cosine  function  as  applied  to  Equations  5.3  and  5.4,  when  the  parameters  of  n',  x, 
and  emissivity  are  fixed.  As  a  result  the  results  of  the  model  differ  from  experimental 
results  in  the  following  ways. 

•  Experimental  parallel  reflectance  values  are  greater  than  those  predicted  by  the 
cluster  model. 

•  Conversely,  experimental  perpendicular  reflectance  values  are  less  than 
corresponding  values  determined  by  the  cluster  model. 

B.  TABULATION  OF  THE  KEY  PARAMETERS 

The  values  for  index  of  refraction  (n'),  cluster  to  bulk  surface  area  ratio  (x),  and 
emissivity  (e)  are  considered  key  elements  in  the  development  of  the  model  since,  it  is 
these  values  which  are  metal  property  depjendent  and  in  large  part  account  for  the  manner 
in  which  light  is  reflected  from  the  metal.  In  this  work  non-spectral  equations  were  used. 
Table  6.1  below,  lists  the  values  of  each  parameter  arrived  at  by  the  iteration  process 
given  in  the  computer  algorithm  development  in  Chapter  V. 
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Table  6.1  KEY  CLUSTER-MODEL  PARAMETERS  OF  METALS 


Metal  Name 

Emissivity  (f) 

Index  of 
refraction 
(n') 

Surface  Area 
Ratio  (x) 

Gold 

0.04674 

2.6696 

0.6074 

Nickel 

0.08800 

3.0710 

0.5663 

Silver 

0.03998 

3.1265 

0.0865 

Aluminum 

0.05216 

5.6479 

0.0607 

C.  ANALYSIS  OF  THE  PARAMETERS 

The  index  of  refraction  for  gold  and  silver,  as  computed  by  the  model  (Table  6.1) 
is  consistent  with  experimentally  obtained  values.  Therefore  this  portion  of  the  model 
development  appears  to  be  validated.  From  Table  6. 1  one  may  note  that  the  value  of  the 
surface  area  ratio  (x)  decreases  as  a  function  of  an  increasing  normal  reflectance.  One 
can  postulate  that  as  the  total  reflectance  at  normal  incidence  increases  this  could  be 
accounted  for  by  a  diminishing  contribution  from  the  cluster,  due  to  less  cluster  surface 
area  thereby  allowing  more  of  the  surface  area  to  act  as  a  metallic  grey  body. 


D.  SIZE  OF  CLUSTERS  VERSUS  PLOTS  OF  "n"  AND  "k" 

Determining  the  size  of  the  clusters  is  an  important  step  in  delineating  which 
portion  of  the  energy  spectrum  would  be  most  susceptible  to  the  cluster  polarization 
model.  Electromagnetic  waves  whose  wavelengths  are  smaller  than  the  size  of  the 
cluster  (i.e.,  X-rays)would  not  see  the  surface  clusters  but  the  interior  bulk  metal. 
Conversely,  electromagnetic  waves  which  are  much  larger  in  wavelength  than  the  cluster 
(i.e.,  long  IR  and  larger)  will  not  be  effected  by  the  cluster  as  the  would  tend  to  be 
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reflected  as  from  an  electric  conductor.  As  postulated  in  Ref.  3,  the  size  of  the  cluster 
is  expected  to  be  around  one  micrometer.  Evidence  of  this  may  be  supported  by 
observing  the  experimental  spectral  plots  of  n  and  k  for  various  metals  over  the  energy 
spectrum.  Tables  6.2  thru  6.5  are  excerpts  from  tables  given  in  Ref.  10  which  will  help 
validate  the  postulate. 


Table  6.2  COPPER  [Ref.  10] 


Wavelength  (/im) 

n 

k 

0.5166 

1.12 

2.60 

0.5390 

1.04 

2.59 

0.5635 

0.826 

2.60 

0.5904 

0.468 

2.81 

0.6880 

0.213 

4.05 

1.240 

0.433 

8.46 

Table  6.3  GOLD  [Ref.  10] 


Wavelength  (/tm) 

n 

k 

0.4592 

1.426 

1.846 

0.5166 

0.608 

2.120 

0.6888 

0.160 

3.80 

0.7749 

0.174 

4.86 

1.033 

0.272 

7.07 

1.240 

0.372 

8.77 

37 


Table  6.4  SILVER  [Ref.  10] 


Wavelength  (/xm) 


0.3155 


0.3237 


0.3306 


0.4769 


1.033 


1.879 


1.044 


0.616 


0.371 


0.132 


0.226 


0.600 


0.514 


0.609 


0.813 


2.72 


6.99 


11.4 


Table  6.5  TUNGSTEN  [Ref.  10] 


Wavelength  (/xm) 

n 

0.03815 

0.797 

0.05209 

0.483 

0.06812 

0.766 

1.265 

3.15 

2.138 

1.18 

6.199 

3.87 

0.286 


0.598 


1.19 


4.36 


8.44 


28.3 


From  the  above  tables  one  will  note  that  the  value  of  n  changes  rapidly  for  all  the  metals 
around  the  0.3  /xm  to  1.0  /xm  range,  except  for  tungsten  which  displays  two  ranges  of 
wavelength  where  n  changes  rapidly.  Further,  from  the  graphs  of  these  plots  [Ref.  10] 
given  in  Appendix  D,  one  will  note  that  while  tungsten  has  two  ranges  of  anomalies  the 
only  anomaly  for  the  other  metals  is  in  the  range  of  wavelength  stated  above.  This 
would  seem  to  indicate  that  some  other  mechanism  of  reflection  is  in  effect  over  the 
range  of  wavelengths  near  one  micrometer,  and  that  mechanism  could  be  a  result  of 
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interaction  of  the  electromagnetic  wave  with  the  cluster  surface.  It  is  stated  in  Ref.  1 
that  in  the  portion  of  the  energy  spectrum  which  is  in  the  near  infrared  xi  is  difficult  to 
measure  the  optical  properties  because  of  the  "anomalous  skin  effect"  [Ref.  9]. 
However,  since  all  the  metals,  except  tungsten,  display  a  similar  variance  of  n  in  the 
same  wavelength  range,  it  is  difficult  to  dismiss  the  possibility  that  the  cluster  size  is 
about  l^m. 
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Vn  SUMMARY  AND  CONCLUSIONS 


A.  MODEL  COMPARISON 

Explaining  polarization  upon  reflection  from  metallic  surfaces  is  not  completely 
understood  [Ref.  1].  As  examples  of  this,  one  need  only  note  that  there  exists  various 
models  each  pertaining  to  a  specific  portion  of  the  energy  spectrum,  which  have  been 
developed  in  an  attempt  to  explain  the  results  of  reflection  rather  than  mechanisms.  For 
instance,  the  Drude  model  (i.e.,  free-electron  theory),  shows  agreement  with 
experimental  results  in  the  long  infrared  and  longer  wavelength  portion  of  the  spectrum 
where  the  interaction  with  the  bound  electrons  can  be  ignored.  However,  this  model 
begins  to  break  down  in  the  near  infrared  and  shorter  wavelength  portions  of  the  energy 
spectrum.  In  fact  Ref.  5  states  just  this  and  refers  the  reader  to  the  quantum  model.  The 
quantum  model  better  describes  polarization  in  the  ultraviolet  region  of  the  spectrum. 
As  stated  in  Ref.  1  this  model  only  partially  explains  the  polarization  mechanism  for  this 
region. 

The  native  cluster  model,  as  formulated,  fills  the  gap  between  the  other  two 
models,  and  expansion  to  encompass  other  regions  of  the  spectrum  seems  to  be  restricted 
by  the  cluster’s  physical  properties  (i.e.,  size  and  the  fact  that  they  reside  on  the  surface). 
Although  still  undergoing  development,  the  formulation  offers  a  mechanism  for 
polarization  in  the  near  visible  and  visible  portions  of  the  spectrum.  From  the  above, 
one  can  postulate  that  the  answer  to  the  fundamental  problem  of  polarization  upon 
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metallic  reflection  will  require  all  three  models  because  of  the  complex  nature  in  which 
different  energy  waves  interact  with  free  and  bound  electrons  of  the  conductor. 

B.  MODEL  EFFECTIVENESS 

Values  of  non-spectral  reflectance,  as  determined  by  the  native  cluster  model,  are 
within  10  %  of  the  measured  values  throughout  the  range  of  angles  of  incidence.  The 
variance  between  the  model  and  experimental  results,  while  not  a  constant,  is  consistent 
in  direction.  That  is,  the  perpendicular  value  of  reflectance  from  Ref.  7  is  somewhat  less 
than  the  value  calculated  by  the  model  (see  Figures  6.1  and  6.2),  and  conversely  the 
parallel  reflectance  from  Ref.  7  is  somewhat  greater  than  that  of  the  model.  This 
consistency  in  the  direction  of  the  variance  gives  rise  to  the  possibility  that  the  model, 
as  it  currently  stands,  is  incomplete.  It  will  be  remembered  that,  from  Chapter  II,  the 
indirect  reflection  from  within  the  cluster  was  not  accounted  for  in  the  calculations. 

However,  the  model  does  generate  curves  with  the  same  characteristics  and  shape. 
More  precisely  the  generated  curves  have  the  same  values  at  normal  and  grazing  angles 
of  incidence  and  the  same  angle  of  incidence  defined  earlier  to  be  the  pseudo-Brewster 
angle.  This  gives  the  model  validity  and  a  measure  of  effectiveness  for  predicting  the 
reflectance  curves  of  "white"  light  from  various  metallic  surfaces. 

The  calculation  of  n'  based  on  the  pseudo-Brewster  angle  appears  to  be  correct  in 
that  it  and  the  index  of  refraction  as  given  by  the  free  electron  model  are  in  close 
agreement.  This  could  become  an  important  issue  since  virtually  all  listed  references 
have  pointed  out  that  surface  preparation  of  the  metal  is  non-standardized  but  yet,  crucial 


41 


in  the  measurement  of  the  optical  properties.  This  fact,  coupled  with  the  statement  from 
Ref.  10  which  says  that  the  measurement  of  the  optical  constants  is  difficult  in  the  range 
above  the  infrared  portion  of  the  energy  spectrum  due  to  the  "anomalous  skin  effect", 
could  make  the  present  model  a  useful  tool  in  helping  to  set  a  new  standard  for  the 
preparation  of  metallic  surfaces. 

The  current  formulation  is  for  non-spectral  visible  waves,  but  could  be  extended 
in  order  to  accommodate  a  spectral  source  of  light.  This  would  be  desirable  from  a 
validation  standpoint,  since  the  author  found  limited  sources  of  data  for  non-spectral 
reflectance  form  metals,  while  finding  an  abundance  of  spectral  data  for  both  reflectance 
curves  and  optical  properties. 

C.  RECOMMENDATIONS 

As  stated  in  Section  A  of  this  chapter,  with  further  validation  and  refinement,  the 
native  cluster  model  could  be  useful  in  filling  the  gap  between  the  free  electron  model 
and  the  quantum  model.  It  is  important  to  note  however,  that  the  development  of  this 
model  does  not  attempt  to  make  other  models  of  polarization  obsolete,  but  rather  offers 
an  alternative  mechanism  for  polarization. 

More  research  needs  to  be  performed  on  bound  native  clusters.  New  work  could 
encompass  development  of  new  observation  techniques.  The  ability  to  detect  clusters  and 
further  observe  their  behavior  over  time  with  varying  surface  energies  is  needed.  This 
would  lead  to  answers  concerning  size,  geometry,  state,  and  perhaps  migration  or  time 
evolution  of  the  clusters. 
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D.  APPLICATIONS 


According  to  Ref.  2,  it  may  be  possible  to  create  made  to  order  clusters  with  the 
most  desirable  electronic,  magnetic  and  optical  properties.  If  the  native  cluster  model 
is  factual,  then  it  may  be  possible  to  alter  the  electron  structure  of  the  clusters  by 
externally  managing  the  energy  balance  of  the  surface  thereby  minimizing  or  maximizing 
polarization  depending  on  the  application.  This  would  be  useful  when  working  with  high 
energy  lasers,  where  a  minimum  loss  of  intensity  and  or  polarization  management  upon 
reflection  from  the  metallic  surface  is  desired. 

Other  possible  applications  would  be  in  thermal  control  design  of  spacecraft  where, 
by  controlling  the  polarization  of  metallic  surface  reflections,  one  could  maximize  or 
minimize  the  absorptance  and  hence  the  emittance  from  the  component  surface.  Using 
similar  methods,  management  of  high  energy  laser  beams  or  enhancement  of  the 
properties  of  off-axis  beam  splitters  would  be  possible.  New  understanding  of  well 
known  phenomena  open  new  vistas  in  applications. 
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APPENDIX  A 


Program  "Reflecta" 
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C  RI^CGQAM  NAME:  REPLECTA 

C  Av'^^CR  SAt-E;  C.  W.  BALOwIH 

C  CE»£-CPE3:  6  APRIL  1991 

C  last  REVISION;  02  DECEHBER  9i 

C 

10  program  REPLECTA 

c 

C  program  REci_ec’A  is  dcsigneo  as  a  toou  POR  the  STUDv 
C  CP  THE  "NATIVE  CLUSTER"  MCOEL  AS  IT  PERTAINS  TO  METALS. 

C  THE  BASIS  CP  this  program  IS  THAT  CLUSTERS  ACT  AS  an  AGENT  IN  The 

C  ^ARTIAu  POLARITA’-ION  OF  VISIBLE  ELECTROMAGNETIC  ;ES  AND  HENCE 

C  SUPACE  AREA  CP  THE  METAL  CAN  BE  OIVIDEO  InTQ  Two  PARTS,  that 

C  CGMPCSEC  GF  The  CLUSTER  AND  TmE  REMAINDER  OF  BULK  MATERIAL. 

C  Sv  combining  the  REF'^ECTANCE  from  both  SURFACE  AREAS  Th£ 

C  total  reflectance  is  calculated  for  parallel  and  perpendicular 

C  CC'-PONEnTS.  the  THEORV  BEHIND  THIS  MODEL  IS  GIVEN  IN  A  THESIS 

C  Bv  C.  W.  BALDWIN  DATED  DEC.  1991 
C 

C  fi£‘'LECT|  IS  THE  DATA  FILE  PROM  WHICH  THE  PLOTS  OF  REFLECTANCE 
C  ARE  generated  BV  INVOCATION  OF  EASV-PlCT,  MOTE  ONE  MUST  RENAME 

C  T-e  reflECTI  file  to  REFuECTl  DATA  IN  ORDER  TO  READ  ThE  DATA  INTO 

c  easv-plot 

c 

C  DECLARE  ALL  VARIABLES 

C 

:0  REAL  PMKIOOJ.  PHlPRlMEdOO).  RSDP(ICO).  RSDSdOOJ 
2'i  real  RBdOO).  RTOTuP  (  100  ) ,  RTOTlS(IOO) 

50  REAL  N.  NPRIHE.  EPSLON.  K 

51  REAL  RTOTLPl(51.Xl(5) 

C 

C  C' EN  THE  FILE  PEFLCTl 

C 

55  OPEN  (UNIT  .  8,  file  •  'REFLECTI-.  STATUS  •  'NEW) 

C 

C  SET  ALL  PARAMETERS 

C 

AO  N  t  1 . 000 

C  N  IS  The  index  of  refraction  of  the  incident  medium  (AIR) 

c 

50  PMiMlNl  1  72.0 

C  PHIHINI  IS  The  angle  of  INCIDENCE  WHICH  GIVES  TM£  SMALLEST 

r  experimental  PERPEnDICUCAR  component  of  REFLECT.-nCE 

c 

52  xs*  0.«667 

c 

c  xs  IS  The  Initial  guess  of  cluster  to  surface  ratio 
c 

55  RMO  •  0.000002200 
C 

0  RHo  IS  the  resistivity  of  the  metal 

c 

5A  STEHF  •  3000.00 

C  STEMP  IS  the  SOV^CE  temperature  (light  bulb,  sun  ETC) 

C 

RSTOTO  •  0.5* 

c 

C  RSTOTO  IS  THE  EXPERIMENTAL  REFLECTANCE  AT  ODEC  OF  INCIDENCE 


c  calculation  of  emissivitv 
C 

C  Rho  »  resistivity  iOHH-CM)  FunTION  OF  TE’-'P .  KELVIN 
C  STE-e  ■  ’■E'*P  OF  illuminating  SOURCE  KELVIN 
C  GCCD  FOR  GREY  OR  BLACK  BODIES 
C 

EPSLON  «  e.5755»((RHO»STEMP)»»0.S) 

C 

C  CalCL'LATICN  of  nprime 
C 

C  NPOIME  [S  THE  theoretical  VALJE  CF  INDEX  OF  HEFBACTION  OF  THE 

c  cluster. 
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PHIMIN  •  (5. 14/18C.0J*  PKlMtMl 
NPRIMC  •  TAN(PNIHIN) 

C 

C  CALCXA7I0N  OF  X  (CLUSTER  TO  SURFACE  RATIO) 

C 

PHlPfllMl  *  ASIN( (M/NPRImE)  •  Q.OSA^J 

RSDPi  •  ((0.0349  -  PMIPRiHl )/(0.0J49*PMlPftIMn 

RBI  .  (  1.0  -  EPSLON  •  0,9994) 

DO  75  J  •  1.2 
DE»Ta  •  0.0S«J 
Xl(J)  *  XS  *  DELTA 

RTOTlPKJ)  •  (1.0  -  XHJ))«RB1  ♦  (Xl(J)ARSDPI) 

75  CONTINUE 
C 

C  CALCULATE  slope  Of  LINE 
C 

R  •  tRTOTLPl(2)  -  RTaTLPl(U)/(Xl(2)  -  XIU)) 

B  •  RTQTlPKC)  -  R  ■  XH2) 

X  •  (RSTOTO  -  Bl/R 

PRINT*.  ‘METAL  TVPE  1$  WICXEL* 

PRINT*.  ‘EPSILON  «  ‘.EPSLON 

PRINT*.  ‘RHQ  a  ‘ .RHO 

PRINT*.  ‘STEMP  »  ‘.STEMP 

PRINT*.  ‘CLUSTER  TO  BULK  RATIO  •  ‘.X 

PRINT*.  ‘NPRIME  •  ‘.NPRIHE 

PRINT* 

C 

C  PRINT  header  for  OUTPUT 
C 

55  PRINT  60 

60  fORNAr(4K. 'PHI'  .5X. ‘PHI  PRIME* .4X. ‘RSDP*  »iX. ‘RSOS* ,5X, ‘RB’ . 

*BX. ‘RTOTLP* .SX.’RTOTLS* ) 

C 

C 

C  INITIATE  LOOP  TO  VARY  PMK 1)  FROM  0-90  DECREES 
C 

70  00  106  t  ■  1.12.2 

72  PMKI)  *  ((-l.O-IJ^CS.U/iaO.O)) 

IF  (PMKI)  ,EG.  0.0)  then 

PHlPHIME(l)*RS0P(I)*RSOS(I)«HBn)*O.OO 

GOTO  1Q4 

EwSe 

CONTINUE 

ENDlF 

•‘'1  PMlPRIMg(I)  •  aSIN  ((N/NPRIME)*(SlM(PHI{n))) 

90  RSDPtn  «  ((TAN(PMHn-PHlPRIME(I)))/tTAWtPMl(I)*PMlPRlHEtni 

•  ))»*2*100.0 

100  KTDsd)  •  ((SlN(PHI(I)-PHlPRIME(n))/(SlM(PHl(I)*PHIPRJMEtn) 
•  ))  ■2"100.0 

101  R9(*)  «(  l.O  -  EPSLON*COS(PHin)))*100.0 

102  RTOT.P(I)  •  (l-Xl-RBd)  *  X*RS0Pd) 

105  RTOT.Sd)  ■  d-XloRBd)  *  X»RSDSd) 

c 

106  CONTINUE 

00  160  I  «  1.92.2 

C  WRITE  data  to  data  file  REFLECTl 

C 

C  CONVERT  PHId)  AND  PHlPRIMEd)  FROM  RADIANS  TO  DECREES 

C 

105  pHtd)  •  pHi(n*dBo.o/s,i6) 

C 

lot  PHlPRIMEd)  ■  PMIPRIMEd)»d«0.0/J.lA) 

RTOTtPd)  •  RT0TLP(3) 

RTOTLSd)  ■  RTOTLSU) 

110  WRITE  (8.120)  PHId  ). PHlPRIMEd). RSDf(I).RSD$tI)>RB(Z) 

♦  .RTOTLPd). RTOTLSd) 

120  fORHATdX.FB. J,IX.F«.S.IX,Fi.3.1X.F».5.IX.F».5.1X.F«.3 

*  .1X.P8.5) 

160  continue 

165  CLOSE  (UNIT  •  8, STATUS  »  ‘KEEP*) 

C 

C  PRINT  OUT  data 
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c 

OPEN  (UNIT  ■  f-PILE  •  ‘ REFLECT  I ‘ . STATUS  •  ’OLD*) 

ISO  DO  1 75  I  •  1 .  <6 

READ  (8.155)  PHltl),  PHIPRIHE(l).  RSOP(l).  R$DS( I ) .RB( 11 . 

•  RTOTLP( I  1 .  RTQTLSt I ) 

155  FOAHAT(lX.F8.3.lK.F8.3.1X.Pa.S.lX.F8.S.LX.F8.S.lX.F8.5 

•  .iX.  F8.S) 

UO  MRlT£(t.l70)  PHKI).  PHIPRIHE(I).  RSOP(l).  RSDSd).  R8(l}. 

•  RTOTLPd).  RTQTLSd) 

170  FORHATdX.F8.S.:x.F8.3.2X.F8.S.2X.F8.5.2X.F8.S.2X.P8.5.2X.F8.S) 

175  continue 
180  END 
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APPENDIX  B 


Reflectance  Data  For  Gold,  Silver,  Aluminum,  Nickel 
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metal  TVPE  is  COLO 

EPSILON  •  0.04(7574 

RHO  •  2.:00000<E-04 

STEMP  ■  SOOO. 0000000 

CLUSTER  TO  BM.K  RATIO  •  0.(075507 

NPRIME  •  2.(45(081 


PHI 

PHI  PRIME 

RSDP 

0.000 

0.000 

0.000 

2.000 

0.745 

20.482 

4.000 

1.457 

20.425 

4.000 

2.244 

20.551 

8.000 

2.588 

20.558 

10.000 

5.750 

20.227 

12.000 

4.447 

20.017 

14.000 

5.155 

15.748 

U.OOO 

5.524 

15.475 

18.000 

4.447 

15.150 

20.000 

7.541 

18.780 

22.000 

8.047 

18.545 

24.000 

8.744 

17.514 

24.000 

5.452 

17.415 

28.000 

10.125 

14.880 

50.000 

10.755 

14.254 

52.000 

11.450 

15.448 

54.000 

12.052 

14.554 

54.000 

12.720 

14.274 

58.000 

15.554 

15.512 

40.000 

15.554 

12.705 

42.000 

14.517 

11.851 

44  '00 

15.084 

10.554 

44.000 

15.455 

10.022 

48.000 

14.145 

5.052 

50.000 

14.478 

8.050 

52.000 

17.171 

7.024 

54.000 

17.445 

5.585 

54.000 

18.055 

4.540 

58.000 

18.525 

5.512 

40.000 

18.555 

2.522 

(2.000 

15.518 

1.555 

44.000 

15.475 

1.184 

(4.000 

20.014 

0.550 

48.000 

20.528 

0.108 

70.000 

20.415 

0.015 

72.000 

20.874 

0.575 

74.000 

21.111 

1.545 

74.000 

21.520 

5.255 

78.000 

21.501 

4.511 

80.000 

21.454 

11.047 

82.000 

21.782 

18.181 

84.000 

21.881 

28.457 

84.000 

21.552 

45.501 

88.000 

21.555 

44.211 

50.000 

22.010 

55.084 

RSDS 

RB 

RTOTLP 

0.000 

0.000 

45.585 

20.720 

55.525 

45.585 

20.777 

55.558 

45.558 

20.871 

55.552 

45.504 

21.005 

55.572 

45.855 

21.178 

55.557 

45.755 

21.550 

55.428 

45.424 

21.(44 

55.445 

45.487 

21.540 

55.507 

45.528 

22.280 

55.555 

45.147 

22.445 

55.408 

48.544 

25.057 

55.444 

48.717 

25.575 

55.750 

48.444 

24.112 

55.755 

48.152 

24.700 

55.875 

47.855 

25.544 

55.552 

47.570 

24.045 

54.054 

47.221 

24.817 

54.124 

44.847 

27.455 

54.218 

44.447 

28.540 

54.514 

44.021 

25.545 

54.415 

45.571 

50.404 

54.524 

45.055 

51.755 

54.457 

44.555 

52.554 

54.752 

44.075 

54.555 

54.871 

45.550 

55.774 

54.554 

42.570 

57.525 

57.121 

42.557 

58.554 

57.251 

41.814 

40.788 

57.585 

41.255 

42.715 

57.521 

40.444 

44.785 

57.441 

40.117 

47.005 

57.804 

55.415 

45.585 

57.545 

55.175 

51.555 

58.057 

58.854 

54.444 

58.247 

58.458 

57.584 

58.555 

58.440 

40.710 

58.555 

58.521 

44.045 

58.705 

55.585 

47.412 

58.844 

40.781 

71.414 

55.025 

42.711 

75.444 

55.185 

45.445 

75.774 

55.544 

50.047 

84.541 

55.508 

54.441 

85.251 

55.470 

45.754 

54.558 

55.855 

75.411 

55.871 

55.554 

55.445 

RTOTLS 
SO. 012 
SO. 012 
SO.OSO 
SO. 115 
SO. 202 
SO. 517 
SO. <58 
SO. 627 
SO. 825 
S1.0<8 
51.505 
51.S85 
SI. 504 
52.2S7 
52.(45 
55.0(5 
55.524 
54.028 
S4.S72 
55.1(2 
55.755 
54.487 
57.228 
58.027 

58.884 
55.810 
(0.802 
(1.8(4 
(5.008 
(4.252 
(5.544 
(4.548 
(8.450 
70.057 
71.775 
75.(05 
75.5(8 
77.(54 

75.885 
82.254 
84.777 
87.455 
50.508 
55.550 
54.552 
55.520 
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hetal  type  is  silver 

EPSILON  •  0.o:y9S2S 

RHO  •  i.*ioeoojE-o» 

STEMP  •  JOOO. 0000000 

CLUSTER  TO  BULK  RATIO  • 

nprihe  •  s.i2«y«s 


PHI 

PHI  PRIME 

RSOP 

0.000 

0.000 

0.000 

2.000 

0.640 

26.536 

4.000 

1.278 

26.473 

6.000 

1.916 

26.370 

8.000 

2.551 

26.224 

10.000 

3.184 

24.037 

12.000 

3.813 

25.806 

14.000 

4.4  38 

25.533 

16.000 

5.058 

25.215 

18.000 

5.672 

24.853 

20.000 

6.280 

24.445 

22.000 

6.882 

23.991 

24.000 

7.475 

23.489 

26.000 

8.060 

22.939 

28.000 

8.636 

22.340 

30.000 

9.203 

21.689 

32.000 

9.759 

20.986 

34.000 

10.304 

20.231 

36.000 

10.837 

19.421 

38.000 

11.358 

18.556 

40.000 

11.865 

17.636 

42.000 

12.359 

16.660 

<4.000 

12.839 

15.627 

<6.000 

13.303 

14.539 

48.000 

13.752 

13.398 

50.000 

14.185 

12.205 

52.000 

14.600 

10. 9« 

54.000 

14.999 

9.686 

56.000 

15.379 

8.375 

58.000 

15.741 

7.046 

60.000 

16.084 

5.715 

62.000 

16.407 

4.410 

64.000 

16.711 

3.162 

66.000 

16.993 

2.021 

68.000 

17.255 

1.049 

70.000 

17.496 

0.336 

72.000 

17.715 

0.006 

74.000 

17.912 

0.234 

76.000 

18.086 

1.265 

78.000 

18.238 

3.453 

80.000 

18.367 

7.311 

82.000 

18.473 

13.603 

84.000 

18.556 

23.485 

86.000 

18.615 

38.757 

88.000 

18.651 

62.307 

90.000 

18.663 

98.953 

0.0864533 


RSDS 

RB 

RTOTLP 

RTOTLS 

0.000 

0.000 

89.998 

90.002 

26.577 

96.004 

89.998 

90.002 

26.639 

96.011 

90.000 

90.014 

26.743 

96.024 

90.002 

90.034 

26.889 

96.041 

90.005 

90.062 

27.077 

96.062 

90.008 

90.098 

27.309 

96.089 

90.013 

90.143 

27.585 

96.120 

90.018 

90.195 

27.907 

96.156 

90.023 

90.256 

28.275 

96.197 

90.029 

90.325 

28.691 

96.243 

90.036 

90.403 

29.156 

96.293 

90.042 

90.488 

29.673 

96.347 

90.048 

90.583 

30.243 

96.406 

90.055 

90.686 

30.869 

96.469 

90.061 

o 

•4# 

CD 

31.552 

96.537 

90.066 

90.919 

32.296 

96.609 

90.071 

91.049 

33.102 

96.685 

90.075 

91.188 

33.975 

96.765 

90.078 

91.336 

34.918 

96.849 

90.080 

91.494 

35.933 

96.936 

90.080 

91.662 

37.025 

97.028 

90.080 

91.840 

38.198 

97.123 

90.077 

92.029 

39.455 

97.221 

90.075 

92.227 

<0.802 

97.323 

90.068 

92.437 

52.242 

97.429 

90.061 

92.658 

<3.782 

97.537 

90.053 

92.890 

<5.426 

97.668 

90.044 

93.134 

<7.179 

97.763 

90.035 

93.389 

49.048 

97.880 

90.027 

93.658 

51.038 

97.999 

90.021 

93.939 

53.155 

98.121 

90.019 

94.234 

55.406 

98.245 

90.025 

94.542 

57.797 

98.372 

90.042 

94.864 

60.336 

98.500 

90.075 

95.201 

63.028 

98.630 

90.132 

95.552 

65.880 

98.762 

90.224 

95.919 

68.900 

98.895 

90.366 

96.302 

72.C9S 

99.030 

90.578 

96.701 

75.471 

99.166 

90.891 

97.117 

79.034 

99.303 

91.350 

97.551 

82.793 

99.441 

92.020 

98.001 

86.752 

99.579 

93.001 

98.470 

90.918 

99.718 

94.448 

98.957 

95.297 

99.857 

96.611 

99.463 

99.892 

99.997 

99.907 

99.988 
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METAL  TYPE  IS  ALUMINUM 

EPSILON  •  0.0S2ISS1 

RHO  ■  2.740000ie-0< 

STEMP  ■  SOOO. 0000000 

CLUSTER  TO  BULK  RATIO  ■  0.0*07 IAS 

NPRIME  •  5.4478*72 


PHI 

PHI  PRIME 

RSOP 

0.000 

0.000 

0.000 

2.000 

0.554 

48.840 

4.000 

0.708 

48.777 

4.000 

1.040 

48.471 

8.000 

1.412 

48.542 

10.000 

1.742 

48.550 

12.000 

2.110 

48.114 

14.000 

2.455 

47.852 

14.000 

2.777 

47.505 

18.000 

5.154 

47.127 

20.000 

5.472 

44.705 

22.000 

5.805 

44.250 

24.000 

4,150 

45.702 

24.000 

4.452 

45.120 

28.000 

4.748 

44.480 

30.000 

5.077 

43.777 

52.000 

5.584 

45.01* 

54.000 

5.482 

42.187 

54.000 

5.774 

41.287 

58.000 

4.257 

40.514 

40.000 

4.554 

57.242 

42.000 

4.805 

38.127 

44.000 

7.044 

54.704 

44.000 

7.518 

55.587 

48.000 

7.542 

54.170 

50.000 

7.77* 

52.448 

52.000 

8.021 

51.015 

54.000 

8.257 

27.257 

54.000 

8.442 

27.580 

58.000 

8.457 

25.570 

40.000 

8.822 

25.224 

*2.000 

8.77* 

20.757 

44.000 

7.157 

18.517 

44.000 

7.511 

15.7*5 

48.000 

7.451 

15.277 

70.000 

7.580 

10.555 

72.000 

7.477 

7.784 

74.000 

7.802 

5.100 

74.000 

7.87* 

2.471 

78.000 

7.777 

0.778 

80.000 

10.044 

0.000 

82.000 

10.102 

1.200 

84.000 

10.144 

4.074 

84.000 

10.178 

17.74* 

88.000 

10.177 

45.574 

70.000 

10.204 

78.187 

RSOS 

RB 

RTOTLP 

0.000 

0.000 

71.777 

48.702 

74.787 

71.777 

48.744 

74.777 

72.004 

47.071 

74.815 

72.012 

47.217 

74.855 

72.024 

47.410 

74.845 

72,057 

47.64S 

74.878 

72.057 

47.720 

74.757 

72.077 

50.240 

74.784 

72.105 

50.405 

75.057 

72.150 

51.014 

75.078 

72.140 

51.447 

75.144 

72.175 

51.747 

75.255 

72.227 

52.514 

75.511 

72.244 

55.110 

75.574 

72.505 

55.752 

75.482 

72.545 

54.445 

75.57* 

72.585 

55.185 

75.475 

72.427 

55.775 

75,777 

72.471 

54.818 

75.887 

72.515 

57.714 

74.005 

72.558 

58.445 

74.125 

72.401 

57.447 

74.247 

72.444 

*0.724 

74.575 

72.484 

*1.845 

74.508 

72.725 

45.017 

74.444 

72,7*0 

*4.250 

74.787 

72.775 

45.545 

74.752 

72.825 

*4.877 

77.081 

72.847 

48.515 

77.254 

72.871 

*7.772 

77.570 

72.88  7 

71.535 

77.547 

72.877 

72.745 

77.711 

72.705 

74.414 

77.87* 

72.705 

74.55* 

78.045 

72.878 

78.142 

78.215 

72.871 

80.054 

78.585 

72.884 

81.778 

78.557 

72.885 

85.787 

78.755 

72.702 

84.045 

78.712 

72.755 

88.211 

77.071 

75.074 

70.424 

77.270 

75.31* 

72.705 

77.451 

75.785 

75.052 

77.452 

74.474 

77.4*5 

77.814 

74.574 

77.7*5 

77.77* 

77.88* 

RTOTLS 
72.001 
72.001 
72.014 
72.055 
72.0*5 
72.104 
72.150 
72.20* 
72.2*7 
72.541 
72.422 
72.511 
72.408 
72.715 
72.827 
72.747 
75.078 
75.217 
75.5*5 
75.517 
75. *78 
75.848 
74.02* 
74.211 
74.404 
74.404 
74.811 
75.02* 
75.247 
75.478 
75.714 
75.757 
74.207 
74.4*4 
*4.72* 
*4.77* 
77.271 
77.552 
77.857 
78.152 
78.450 
78.755 
77.041 
77.554 
77.471 
77.775 
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HETAL  TYPE 

IS  NICKEL 

EPSILON  • 

0.0447S74 

RMO  • 

2.2000004E-04 

STEMP  • 

SOOO. 0000000 

CLUSTER  TO 

BULK  RATIO 

0. 

.5442929 

NPRIME  • 

S.0710: 

201 

PHI 

PHI  PRIHE 

RSDP 

RSDS 

RB 

RTOTLP 

RTOTLS 

0.000 

0.000 

0.000 

0.000 

0.000 

55.989 

54.012 

2.000 

0.451 

25.859 

25.900 

9S.S29 

55.989 

54.012 

4.000 

1.S02 

25.798 

25.942 

95.SS8 

55.958 

54.051 

4.000 

1.951 

25.495 

24.045 

95.S52 

55.904 

54.115 

•  .000 

2.597 

25.551 

24.210 

9S.S72 

5S.8S2 

54.204 

10.000 

S.241 
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